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I calculate the effect of interactions among superconducting quasiparticles in two-dimensional(2D) 
a superconductor at T=0. Tlie strength of the effective interaction among the quasiparticles is es- 
sentially given by the screened Coulomb interaction which has strength at low frequency because 
^ ' of the gapless nature of the plasmon. This is in contrast to three dimensions where the effective 

interaction has negligible weight at frequencies ~ A, the superconducting gap. The quasiparticle 
interactions give rise to strong-coupling effects in experimental quantities which are beyond the con- 
QQ ' ventional Eliashberg treatment of superconductivity. The present calculation offers an explanation 

of why these effects are much larger in 2D than in 3D superconductors and, in particular, why the 
analogous strong-coupling effects due to quasiparticle interactions are seen in data on the quasi-2D 
cuprate superconductors, the strong-coupling features seen in data on the cuprates are discussed in 
light of the present calculation. 
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The BCS theory of superconductivity and its extension Eliashberg theory describe the electronic properties of the 

superconducting state in terms of self-energies which arise from the interaction of the superconducting quasiparticles 

Q 1 with a boson which provides the mechanism for superconductivity [0,0. It is assumed that the direct Coulomb 

^ interaction among the quasiparticles leads to renormalizations of the quasiparticle effective mass and lifetime in 

the normal state but that these renormalizations are unaffected by the superconducting transition. Within this 

approximation direct interactions among the superconducting quasiparticles are ignored and the only indication 

^ , of interactions in the superconducting state are the a'^F{Lo) strong-coupling features which appear in tunneling 

00 ■ experiments and measurements of optical conductivity. The a^F(uj) function then gives information on the density of 

Q^ , the bosons and the strength of interaction with the quasiparticles. This has proved to be a very successful treatment of 

~^ the properties of conventional superconductors where this mean field approximation for the direct long-range Coulomb 

I . interaction is justified by data. Models for the mechanism of high Tc superconductivity are based on the magnetic 

[~^ ■ properties arising form strong short range Coulomb correlations as described in the Hubbard model. Treatments of 

On ■ these models also frequently assume that the normalizations determined in the normal state are unaffected by the 

transition to superconductivity. 

However quasiparticle interactions in the ordered state have been shown to renormalize the superconducting gap, 
change the temperature dependence of the gap, and lead to finite quasiparticle lifetimes &-0|. In the cuprate super- 
conductors quasiparticle interactions lead to a feature in tunneling conductance across an superconductor-insulator- 
P^ ■ superconductor((7s/5) junction at 3A in addition to the peak at 2A expected from the mean field approximation [^,||. 
O ■ It is the scaling of the dip feature with A which distinguishes these corrections from the usual features in a^F{Ld) 
^ ' associated with the electron-phonon interaction. This scaling was originally demonstrated by comparing the tunnel- 
ing conductance curves for a set of cuprate superconductors whose Tc's varied from 5K to lOOK [gj. More recently 
the intrinsic nature of this strong-coupling feature was demonstrated by comparing the tunneling conductance for 
Bi2Sr2CaCu208+5 as a function of doping and seeing the dip feature scale with the position of the 2 A peak ^. The 
C^ ' size of the dip feature is much larger than the conventional electron-phonon corrections in the cuprates. This dip fea- 
ture has been seen in the tunneling results of a number of groups ||lQ-n3 and the corresponding feature in the spectral 
density has also been identified in ARPES(Angle Resolved Photoemission Spectroscopy) data on Bi2Sr2CaCu208+5 
[p"5|-p7[ and has recently generated a renewed interest ||l^,|l^ • 

The magnitude of these corrections to the mean field treatment of superconductivity in the cuprates distinguishes 
the cuprates from conventional superconductors. Corrections to the mean field approximation, although present 
in conventional superconductors, do not in general lead to experimentally detected features. In this paper I argue 
that it is the low-dimensionality of the cuprate superconductors which is responsible for the comparatively large 
effects of quasiparticle interactions in the cuprates compared to conventional 3D superconductors. In low dimensional 
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systems fluctuations lead to the destruction of long-range order. In one dimension quantum fluctuations of the order 
parameter are sufficient to destroy long-range order pC| and in two-dimensions thermal fluctuations destroy the order 
pl| , p2[ . Given this the effects of quantum fluctuations on superconducting properties would be expected to be more 
important in 2D than in 3D. This is investigated using a model Hamiltonian for 2D and 3D s-wave superconductors 
at zero temperature which includes the Coulomb interaction between electrons and the pairing interaction originally 
introduced by BCS which is characterized by a magnitude Vq for electrons within an energy cod of the Fermi surface. 

In this paper I show how the long-range nature of the Coulomb interaction leads to a qualitative difference in 
the quasiparticle properties of 2D and 3D s-wave superconductors. As Anderson ||2^ pointed out the long-range 
nature of the Coulomb interaction leads to the disappearance of the gapless mode associated with the phase of the 
superconductor order parameter in three dimensions and collective mode becomes the plasmon at energies much 
higher the energy scale of superconductivity. As a result the collective mode has no influence on superconducting 
properties p4|~p7[. This is not the case in two dimensions where the collective mode remains gapless just as the 
plasmon mode does for the 2D electron gas ||2q-|3^. This leads to the possibility for the collective mode having a 
strong influence on quasiparticle properties unlike the case in three dimensions. Interactions among electrons are 
described by the Coulomb interaction and the pairing interaction responsible for the s-wave superconductivity in the 
generalized random phase approximation. I show that the quasiparticle interactions in the 2D superconductor are 
stronger than in the 3D case and that this difference can be traced to the low lying collective mode and the enhanced 
magnitude of the effective interaction at low frequencies compared to that in 3D. 

Most models of the electronic properties of the cuprate superconductors have emphasized their layered nature and 
have employed 2D or quasi-2D Hamiltonians. These models have stressed short-range correlations either through an 
on-site Hubbard repulsion, in the weak or strong-coupling limits, or phenomenologically through residual short-range 
antiferromagnetic order which survives in the doped materials |3qJ3q] . In these models a tight- binding bandstructure 
is used with either an on-site Hubbard repulsion coming from the weak-coupling side or the t — J model coming from 
the strong-coupling side. Both approaches have been shown to be consistent with the (1^2 _y2 symmetry of the order 
parameter. The importance of long-range correlations coming from the Coulomb interaction in quasi-2D systems, 
such as the layered cuprate superconductors, have not been considered in these treatments but potentially have a 
strong influence on the properties of quasi-2D superconductors. 

The single-particle self-energies are calculated here for 2D and 3D s-wave superconductors and it is found that their 
magnitude in 2D is roughly ten times the value in 3D. The single-particle self-energies also have a stronger frequency 
dependence in 2D than in 3D which leads to stronger signatures of quasiparticle interactions in the spectral density 
and in the tunneling conductance of 2D superconductors compared to those in 3D superconductors. I calculate the 
strong-coupling corrections in model s-wave superconductor introduced and discuss how these results are related to 
experimental data on the cuprates. The paper is organized as follows. In section II the model is introduced. In section 
HI the nature of the collective modes and their coupling to quasiparticle excitations are reviewed. The calculated 
single-particle self-energies are discussed in section IV and the influence of quasiparticle interactions on measured 
quantities is discussed in section V. The implications of the results presented here for work on the properties of the 
cuprates are discussed in section VI. A preliminary report of this work has appeared elsewhere |37l| . 

II. MODEL 

The starting point of the present analysis is a Hamiltonian, Eq.(l), describing a system of fermions interacting via 
a potential, U{q). 

H = y^£rci cr -f- V V(k,q]ci cK c .-, c.- + - V U(q]ci cK c r.aCt (1) 

A^^^k ka ka 2 ^-^ ^ ' k-qa ~kl+q-a -kl-a ka 2 ^-^ ^ ' k-qa -kl+qP "kip ka V ; 

k,a kf,k,q,a kf ,k^q^a,(3 

where £r = U t^, kp is the Fermi momentum, Uiq) = ^n^, Cs is the static dielectric constant of the material 

and V{k,q) is the same pairing interaction as introduced by BCS, V{k,q) = — VoOlt^D — I^fel)0(w'i3 — ICfe_fl); with 
Q{x) — for a; < and = 1 for a; > 0. wd is a cutoff energy for the pairing interaction, which is a free parameter in 
the current calculation, was taken to be the Debye energy by BCS. Assuming that the ground state of the system at 
low temperatures is an s-wave superconductor, there is a transformation of the Hamiltonian to one written in terms 
of the quasiparticle operators, 7^:^, which destroy this state. Transforming the Hamiltonian it is written in terms of 
these quasiparticle operators, 7I and 7^^, the quasiparticle spectrum, Ej; — y^^p + A^ , and the coherence factors. 



u?. = 1 — Tj?. = i(l + ■^). A is determined by the weak coupling gap equation. The Hamiltonian becomes 
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The terms in the Hamihonian from the Coulomb interaction are, 



Hc"^^ = J2 Uiq)miki + q, fci)m(fc2 + q, h^t^ 



+q]^-kiC-'^2i^k2+q]' 



ki,k2,q 
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^c ' = Z^ ^/(^^^(A:! + (f, fci)m(fc2 - 9, ^2)7l^_-^7£^-^7L^:^^7fe,, + h.c. 
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4'^-^ E ^('?^™(fc^i+'?''^iMfe-9lfc2)7,\^^-^7^,7t^.^^7l,-^,+/^-c. 
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In the 7/^ the net momentum of the two quasiparticles is the same as that in the Coulomb interaction. The H^ 
term describes the scattering of pairs of quasiparticles whose net momentum Q is different from that in the Coulomb 
interaction and appears in the calculation of ladder diagrams with the Coulomb interaction. The Hq term describes 
the creation or destruction of two pairs of quasiparticles one with momentum g and the other with momentum —(fin 
which g is also the momentum carried in the Coulomb interaction. The functions l{ki,k2), m{ki,k2), n(ki,k2) and 
p{ki,k2) are the usual combinations of coherence factors M, 



m{k + q,k)= u^^^v_f: + vj:^^u_j:, 
n{k + q,k)^ uj:^quj: - v^^^v^, 
l{k + q,k)= u^+^u_^ + Vj^^JV_^, 
p{k + q,k)= uj:^^v_j: - v^^M_^. 

The terms from the pairing interaction are, 
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^ki+q]^-kii ^ ^-k2i"^k2+q^ 
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The Hp "'"^ and H^"^^' terms generate phase and amplitude fluctuations of the order parameter. 



In normal ordering the 7I and 7^^ operators terms quadratic in the new operators as well as terms without operators 

are generated from the kinetic energy and interaction terms. Considering first these terms from the kinetic energy and 

pairing interaction they include the effect of the pairing interaction calculated in the Hartree approximation. Ignoring 

the shift in the chemical potential due to the pairing interaction and using the weak-coupling gap equation these terms 

give the kinetic energy term, ^r Er'^l 7?: , in the transformed Hamiltonian. Consequently in calculating the effect 

of the interactions in Equation(2) all contributions which include Hartree contributions to the quasiparticle lines from 

the pairing interaction should be dropped to avoid double-counting. 

The quadratic terms arising in the normal ordering of the Coulomb parts of the Hamiltonian lead to the Hartree- 

2 7, 
Fock spectrum. In the normal state in 3D this leads to a large shift in the chemical potential, (5/i = — "^ ^ , and a 

logarithmically diverging effective mass |3^]. In 2D the self-energy leading to the Hartree- Fock spectrum is E^(fc) = 
-'^''^'"' E{-^) ior k < kp imdY^^D^k) = =^^^[E{^)-{l-{^f)K{^)]iork> fc^, where ^(a;) and £;(x) are elliptic 

integrals of the first and second kind. As in 3D there is a large shift in the chemical potential, 5^ — — I±Ji£. and the 
effective mass diverges logarithmically. These results for the Coulomb interaction in the Hartree Fock approximation 
are strongly effected by higher order terms in the perturbation expansion and the long range interactions should be 
screened Q,Q . The statically screened Coulomb interaction leads to a smaller shift in the chemical potential and to a 
modest reduction in the effective mass. In the calculations discussed in the present paper the normal state quasiparticle 
spectrum is measured from the renormalized chemical potential and the reduction in the effective mass due to the 
statically screened Coulomb interaction is ignored. The single-particle self-energies for the 7I 7^ and for the 7I 7 - 

propagators are calculated at zero temperature using the generalized random phase approximation(GRPA). In this 
approximation the effective interaction is essentially the screened Coulomb interaction minus its static limit. 

III. THE EFFECTIVE INTERACTION AND COLLECTIVE MODES 
A. The Effective Interaction 

The single-particle self-energies are calculated at zero temperature using the generalized random phase approxi- 
mation(GRPA). In the GRPA ring diagrams, where each ring is the two quasiparticle propagator, are summed in 
which the vertices are dressed by the pairing interaction. In the absence of vertex corrections the two-quasiparticle 
propagator is 

r q — kk q — kk ^ q — k k' 

The interaction connecting the dressed two-quasiparticle propagators carries the net momentum and frequency of 
the two quasiparticles and is described by terms Hq and Hp in the Hamiltonian. However the magnitude of the 
pairing interaction is negligible compared to the Coulomb interaction and so Hp can be ignored. Scaling energies 
with respect to A and momenta with respect to pp, the strength of the interaction is given by a dimensionless 
parameter ^^!^-7V(0) — — - — , where as is the Bohr radius. The strength of the effective interaction is determined 
by the product tspp which is a free parameter in the model. 

At zero temperature the vertex corrections in the rings are given by the repeated scattering the two quasiparticles 
which is described by the Hp °"'^'^' and the Hp"^^' terms in the Hamiltonian. Calculating a "rung" in this ladder sum 
by summing over the intermediate quasiparticle momenta, one finds 

B2o{q,uj)l{q~ ki,ki)l{q- fc2,fc2) + 6*20(9, w)n((f- ki,ki)n{q- k2,k2) (13) 



+ £^20(9,^) 
where 
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1 , ^j: , C--fc , ^n-k + ^k 



\~^ J- , Sife I ^q—k ^ q — k 

fc 



^-(^-)-1.2(^'^_..^._(^_. + ^^). 



Summing over /c, E2o{q, lo) vanishes with particle-hole symmetry. i?2o('fi '^) and C2o('f, i^) are the phase and amplitude 
fluctuations of the order parameter, respectively. As a result of particle-hole symmetry these two channels are 
decoupled from one another and the ladder sum becomes 



-/((f- fci, ki)l(q- fc2, fc2) -I- .. „ ,^ — ^n{q- ki,ki)n{q- fc2, ^2) 



\+fciT^fca'^9-fc2T'^fc2i- 

(15) 



1-^0^20(9,^) ' ' 1 - VoC2o{q, i^) 

This decoupling of amplitude fluctuations from charge fluctuations has been noted previously by Wu and Griffin M . 
The two different combinations of coherence factors in the phase and amplitude fluctuation channels, l{q— fci, fci) 
and n{q — fci, /ci), lead to different couplings to the H)j vertex which are given by a function 020(5, w), for the phase 
fluctuations, and by d2o{q,Lo), for the amplitude fluctuations, where 

f^ \ STt ^^ \ ^q-k + ^k ,.„. 



0^20(9,^) 






The coupling of the H^ to amplitude fluctuations of the order parameter is zero because of particle-hole symmetry. 
As discussed above the static part of the effective interaction leads to a renormalization of the chemical potential and 
the Hartree-Fock spectrum and is to be dropped assuming that energies are to be measured from the renormalized 
chemical potential. The effective interaction becomes 

'^^^^'^'"^ 1 - Umi<l,u:) 1 - U{q)Piq,u = 0) (^'^ 

where 

n(,-,c.) = Ao„(g»--i|i2l^ (18) 

1 - VoB2Q{q,uj) 

Since the weak coupling gap equation, which determines the magnitude of A in the present calculation, is 1 = 
B2o{q~ 0,uj — 0), l — VoB2o{q,i^) becomes VoSB2o{q,i^), the Vo factor cancels in Eq. (18), and there is no the explicit 
dependence of Il{q,uj) on Vq. Furthermore since the energies of interest here are < lOA, esPF is the most important 
parameter and A sets the energy scale. 

B. Collective Modes 

Within this approximation, the imaginary part of the effective interaction is 

V:'fjiq,u;) = U^{q)x\q,oj) = |^ ^'J(]"i^^|.^^)|2 + U{m^ U{q)IiA{q.u:)] (19) 

The second term gives the collective mode contribution. Collective modes have been investigated extensively for 
3D superconductors p3|-p5[ and more recently for 2D and quasi-2D superconductors ]29|-p3[. As is well-known the 
energy of the collective mode in 3D superconductors is comparable to the plasma energy and the collective mode 



has no influence on properties at energies ~ A. On the other hand lj^ = \/ ^'^a^ 9 ^^ long wavelengths in 2D 
superconductors, where N(0) — ATrm/h'^. As a result Veff{q,oj) has an imaginary part given by nZ^biu — lo^) for 
values of c^ - A, where % = [[/2(g) 9nA(^]-i|^^^^^^ 

Vg'f f ((f, w) is U'^{q) times the screened response function which is the dynamic structure factor, S'(g, w), calculated in 

the random phase approximation. Collective modes exhaust the /—sum rule, / duJU!S{q, uj) = 3^, at longwavelengths 
so that 

Nq"^ 
dujujS{q,uj) ~ Z„uj„ = — — (20) 

2m 

2 -. 
N is the number density of particles and the weight in the collective mode is Zq ~ A^^^ in the long wavelength 

limit. As Lo^ — )■ 2A, lo^ increases less rapidly than ^ and is Landau damped in the two quasiparticle continuum. 
This is shown in figure 5 for esPF — SttA^^. At the same time Zg- also starts to increase less rapidly than q^ and 
eventually goes to zero as uj^ reaches 2A. 

In Figure 1 fq is the fraction of the /—sum given by the collective, fq — Z^/^^. When w^- — > 2 A, fq — > 0. At 
larger wavenumbers, the collective excitation reemerges as a resonance above the continuum. In the normal state 
when the collective mode, the plasmon, is calculated in the random phase approximation, the excitation is a pole of 

the response function at the same wavenumbers. However in the present calculation the imaginary part of AoQ^q^co) 

* 2 4. 

behaves as ~ A^(O)-^ at high frequencies. This is because the 7^:^ operators are mixtures of c^:^ and cl and so 

v4oo(g, w) is a mixture of both normal state particle-hole and particle-particle propagators. Consequently there is a 
resonance at the plasmon energy rather than a pole. 

In Figure 2 I show the plasmon resonance in x {q^^) s-s a function of (f for espp — SttA^^. The important point 
about the effective interaction in two-dimensions is that there is weight at low frequencies comparable to A at long 
wavelengths. For example with espp — SttA"^, oj^ < lOA for \q] < 0.125p_F. As a result interactions between 
quasiparticles at low energies ~ A are much stronger than the interactions between the corresponding quasiparticles 
in 3D superconductors. As the value of esPF decreases the strength of the interaction increases and the resonance will 
emerge from the continuum at smaller q and uj. I will discuss the dependence of the results on esPF below. 
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FIG. 1. Collective mode and weight in f-sum for A = 0.0273ei?. and e^pp — SnA ^. 
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FIG. 2. Collective mode emerging from the continuum seen as a resonance in the imaginary part of 
the response function, x il,^)- As q increases the resonance appears at higher energies compared to A 
and becomes more like the normal state plasmon. Eventually it is Landau-damped in the quasiparticle 
continuum just the normal state plasmon. 

Fertig and DasSarma p9| have calculated collective modes in layered superconductors allowing for hopping between 
planes. They found that collective modes exist in the gap for a limited range of momenta in the layers, |<7|, and 
transverse to the layers, qz- This region of \q\ and qz values shrinks with increasing strength of the interlayer hopping, 
W. Interlayer hopping also leads to a gap in the collective mode dispersion at long wavelengths comparable in 
magnitude to W. The important feature of the collective mode spectrum from the point of view of the present work is 
that they are at energies of order A. As long as this is the case low energy quasiparticle interactions will be enhanced 
above those in 3D superconductors. Estimates given by Fertig and DasSarma suggest that this is the case in the 
cuprates. A detailed account of quasiparticle interactions in layered superconductors is in preparation p2[. 



IV. SINGLE-PARTICLE SELF-ENERGIES 



The single-particle self-energy for the 7I 7^ ^ and 7_j: i 7j: -^ propagators are calculated at zero temperature 
using the generalized random phase approximation(GRPA). The self-energies associated with the two single-particle 
propagators, G^t^{p,iEn) and G.y^{p,iEn) are. 



T.^i^{p,iEn) = -7^Yl 



, si/fi^Q 



tanh{——)[n (p - q,p)Ve{q, iE^ - %-g) - m {p- q,p}Ve{q, -iE,, - £^-j?-g)] (21) 



2T 
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lEn — LO — Ep^ql lEn + U) + E^f^g 



tarLh{^^)m{p,p-q)n(p,p-q) x 



q 

/°° duj ,/ w ,,.,,, , 
— coth{ — )V:\q,u;) 
.00 TT 2r E^^ + {u + Ep^g^y 



U! + Eff-g- 



Veiq, iEn - EjJ-g') - Veiq, ~iEn - £^p-g) 



(22) 



Analytically continuing to the real axis and taking the zero temperature limit, the real part of the retarded self- 
energies are 



s;t>-'.^) = -|E 



n^{p- q.P) Vi{q, -E - E.p_^) 



I nOQ 1 



du; 2{E-Ep^,-)V:'{q,iu) 
TT (E- Ep_^Y - co2 

duj2{E + E_j;_^)Vl'{q,uo) 



[E + E^p^,-) 



(23) 



^'■y-yiP^E) = --^m{p,p-q)n{p,p-q) 

X K' (9, E - Ej,^^) - K' (g, -E- Ef!^^) + / — 



2{E - Ep_^)Vi'{q, Lo) ^ 2{E + E_p_^)Vi'{q, u; 
{E - Ej;_^f - u^ 



{E + E^p^^f 



The imaginary parts of the self-energies are, 



p— q^p 



(24) 



(25) 



where V^.^ Aq,Lo) and Vg"irf((f, w) are the real and imaginary parts of the effective interaction. 

The calculated imaginary part of T,^',^{p, E) is shown for a two-dimensional s-wave superconductor in Figure 3 with 
tsPF equal to SttA"^, 
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FIG. 3. Imaginary part of T,^f^{p,E) for a 2D s-wave superconductor for different values of Ep with 
£sPF = 25.13A~'^, LUD = O.leF and A = 0.0273eF. The dashed line is the contribution from the collective 
mode for each Ep. This is very small for Ep = A (dot-dash curve). 



In Figure 3 the dashed curve is the contribution due to scattering off the collective mode and the solid line is the 
sum of the collective mode and continuum contributions. In 3D S''^ (p, E) increases monotonically with increasing 
E > 3Ao as the phase for decay into quasiparticle pairs increases. Comparing the 2D and 3D results one notes that 
the collective mode contribution is absent from the 3D case and that S''^ {p, E) is an order of magnitude larger in 
2D. There is little dependence on \p\ for \p\ > l.Olpp in the continuum contribution to S''^ {p, E). This contribution 
is a monotonically increasing function of energy at these low energies. The different coherence factors which appear 
in the expression for S" (p, E) lead a smaller magnitude than for Y," ^ {p, E). However the energy dependence in the 
collective and continuum contributions arc unchanged. 




FIG. 4. Real part of the self-energy, Ti^^^(p, E), for a 2D s-wave superconductor for different momenta 
with esPF = 25.13A"\ ljd = O.leF and A = 0.0273eF. 



The real part of the Y^^i^ip, E) self-energy is plotted for esPF — SttA"^ in Figure 4 and for €sPf = 47rA~^ in Figure 
5. From Figure 4 one see that the E dependence is different for E < Ep where the dependence is non-monotonic, for 
E "^ Ep where there is very little dependence on E and for E > Ep, after a rapid increase, there is monotonic almost 
linear increase with E. For E ~ 2A this rapid increase has a step-like dependence. 

In Figure 5 Ti^t^^p, E) is seen to have the same form for different values of SsPf but that the magnitude increases as 
EsPf decreases, that is as the magnitude of the Coulomb interaction increases. For Ep > 4A the on-shell self-energy 
'E^f^{p,E — Ep) is only very weakly dependent on E and the magnitude is larger for smaller values of e^pF- For 



^7 ' 7 

esPF — Sir A" 
such that S' 






Ep for small E.^ 



0.3A,and for esPF = AttA'^ E' ^ (p, Ep) -^ -0.78A, as -B -> lOA. For values of esPF 
and corrections to GRPA should be calculated. A detailed discussion of the 



influence of the form of Yi^t^^p, E) on the density of states and the tunneling conductance is given in the next section. 
The other self-energy, S' (p, E), is shown for 2D in figures 6. The magnitude of this self-energy is smaller than the 
magnitude of E' ^ {p, Ep) ~ Ep and falls off monotonically for large E. Its magnitude also increases with decreasing 
esPF but its dependence on E remains unchanged. 

Comparing the 2D and 3D cases again one finds that the magnitude of the real part of the self-energies is roughly 
an order of magnitude smaller in 3D than in 2D and for E > Aq the energy dependence is almost linear without the 
step feature of the 2D case. The same free parameter, esPF, is present in the 3D calculations however because of the 
magnitude of the self-energies in 3D -^ remains small for much smaller values of CsPf- 




O.O 1 -O 2.0 3.0 4.0 5.0 e.O "7.0 SO 



FIG. 5. Real part of the self-energy, T,^f^{p, E), for a 2D s-wave superconductor for different momenta 
with lod ~ O.leF and A = 0.0273eF- for eapp = A-kA"^. The dependence on E is the same as that for 
esPF = SttA^^ but the magnitude is three times larger. 



0.100 



0.050 




FIG. 6. Real part of the self-energy, Ti^~,(p,E), for a 2D s-wave superconductor for different momenta 
with tsPF = 25.13A"\ LOD = O.lcF and A = 0.0273eF. 
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The same calculation can be carried out for a 2D d-wave superconductor. In many models of superconductivity 
of the cuprates the d-wave symmetry arises from a pairing interaction with strong short-range correlations which 
are not present in this work. In a more complete treatment of the charge carriers in the cuprates both long range 
and short range parts of the Coulomb interaction would be included. As Hwang and Das Sarma |2|] have pointed 
out the long range Coulomb interaction leads to quantitatively the same collective mode behavior in both s- and 
d- wave superconductors so that there is the same kind of enhancement of quasiparticle interactions in a 2D d- 
wave superconductor as appears in the 2D s-wave case. In calculations of superconducting properties for a d-wave 
superconductor with a tight-binding bandstructure bandstructure and gap anisotropy lead to anisotropy in the self- 
energy Q. Nodes in the gap function also lead to damping of collective modes. These effects will not lead to a 
qualitatively different conclusion from the simple model discussed here but are clearly important for comparison with 
data on the cuprates. 



V. STRONG COUPLING FEATURES IN TUNNELING CONDUCTANCES 



Considering the expression for the the tunneling conductance across an superconductor-insulator-superconductor 
junction, 

u;)\T\^ (26) 



g 



dI{V) 



dV 



-y 



dV 



dujA{ki,uj)A{k2,eV - 



fclfc2 



where I(V) is the current across an SIS junction and T is an unspecified matrix element between states on opposite 
sides of the junction which in principle depends on ki and fc2. In this calculation the tunneling matrix elements are 
taken to be constants Es]. In describing tunneling in systems with strong bandstructure effects this approximation 
for the tunneling matrix element is not always justified p4| . Summing the carrier spins the spectral density A{k,uj) 
is 



A{k,E) = - 



ulKiAk.E) 



k in 



vlY.'' (- 



^E) 



[E- 



Ek 



S'tJfc,i?))^ 



(s;'t^(fc,i?))2 



[E- 



^k- 






^E)Y + {i:':^,^{-k,-E)Y\ (27) 




FIG. 7. Comparison between A'^(e) with the mean field approximation, A'^ (e) for 2D superconductors. 
jyAiF^^-j was calculated using a Lorentzian with To — A/20. 
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FIG. 8. Comparison between gsis corresponding to the N{u>) in Figure 7(a) with the mean field approx- 
imation gg'j^ (dashed curves), for 2D superconductors. N^''^(e) was calculated using a Lorentzian with 
Fo = A/20. 

The tunneling density of states, essentially the tunneling conductance across a supcrconductor-insulator-normal 
nietal(SIN) junction, and the tunneling conductance are compared with the mean field approximation in Figure 7 
and in Figure 8 for esPr = &TrA~^. There is a piling up of states at E's ~ A due to T,' ^ (p^E) compared to the 
density of states calculated in the mean field approximation, N^^^{uj) == "^^ , _p°)2_^-n'i - In these calculations an 

extra damping term, Fq = ^, is added to S"^ {Pi E) to mimic the effect of impurity scattering on these features. 
This piling up of states leads to a non-monotonic dependence of the current across an S-FS junction as a function of 
applied voltage which in turn results the tunneling conductance shown in Figure 7. The large variation in gsisi^^) 
at eV ~ 2 A comes from the magnitude of E' ^ {p, E). 

The detailed E-dependence leads to strong-coupling effects which are seen at different frequencies in N{uj) and 
have corresponding features at frequencies shifted by A in gsisisV). These strong-coupling features are shown in 
the density of states and in the tunneling conductance in Figures 9 and 10 for different values of a phenomenological 
damping factor Fq in the propagators. In this way some idea of the sensitivity of these features with damping can be 
seen. The magnitude of these features are at the 5% level in the density of states for Fq = A/20. The most prominent 
of these features are a bump feature between ~ 2.6A and 2.9A and the step feature at w = 3A in N{lu). These 
lead to more pronounced features at eV ~ 4 A in gsisi^V). This step feature is the analog of the dip feature seen 
in gsisisV) at eV = 3 A in the d-wave cuprate superconductors whereas the bump feature comes from the collective 
mode contribution as will be discussed below. The origin of these features is the self-energy arising from quasiparticle 
interactions mediated by the screened Coulomb interaction and not from the pairing interaction. 

The on-shell imaginary part of the self-energy, S"| (p, Ep), is given to a good approximation for e^pp ~ SttA^^ by 

S"t^(p, Ep) = Fo + 0.055e(Sp - 3A)[Ep - 3A] (28) 

The small magnitude means that the spectral density at any momentum is to a very good approximation a Forcntzian 
and there is no structure from the energy dependence of 'E^i^{p, E) at values of E where the spectral density, A{p^ E), 
has any appreciable magnitude. If S"^ (p, E) is set equal to Fq = ^ and the conductance recalculated the same 
features appear in the tunneling conductance at the same energies. The only difference is that the new curve lies below 
the curve calculated with S''^ (p, E) . So the E dependence of Y/\ (p, E) is not responsible for the strong-coupling 
features in the density of states. 
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FIG. 9. Strong Coupling effects in the Density of States shown in Figure 7. The sohd curve is the 
density of states calculated with Fo — A/20, the heavy dashed curve is calculated with Fo — A/10 and 
the dotted curve is the mean field density of states with Fo — A/20. 




FIG. 10. Strong-coupling effects at oj = 4A in gsis corresponding to that a,t uj = 3A in the density of 
states shown in Figure 7. The solid curve is the case of Fo — A/20, the dashed curve is for Fo = A/10 
and the dotted curve is the mean field curve with Fq — A/20. 
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The real part, S' ^ (p, E), can be considered to the sum of two contributions, the on-shell piece, S' ^ (p, Ep), and 



I'l 



the energy dependent part, ST,' ^ (p, E). The position of the peak in the spectral density is given by 



7'7' 



Ep + j:'Jp,Ep) + SJ:',Jp,u;p) 



(29) 



ujp is the rcnornializcd quasiparticle spectrum. S' ^ {p,Ep)/A is plotted versus Ep/A in Figure ll(curvc A). 
S' I {p, Ep) has no strong dependence on Ep except in the region near Ep = 3A where the the features are seen 
in the density of states. Beyond Ep — 2A it is a slowly increasing function of Ep except for these fluctuations. 
t^p — Ep is plotted against ujp in the same figure(curve B) and it is seen to have the same rapid variation near 3A. 
Comparing these curves with the density of states in Figure 9 the onset of the bump feature in the density of states 
at u; ~ 2.6A and the step feature at 3A can be seen to arise from structure of the quasiparticle spectrum. This 
structure in cjp causes fluctuations in the density of states of the renormalized quasiparticle spectrum. Consequently 
the strong-coupling features in the density of states and the SIS tunneling arise from the momentum and energy 
dependence of S' ^ {p, E) due to quasiparticle interactions. 

The form of S' ^ [p, ujp) conies from the continuum and collective contributions to the imaginary part of the sclf- 

^cont I 

self-energy is 



7'7 



energy, S'^o"* (j)^ E) and S^° {P^ E) ■ From the Kramers-Kronig relation, the energy dependent real part of the 



S't^(p,^p) = 



1 /-^-S-"* {p,E')dE' 



"7^ 



^'^'T'-^' ^^ TT./o UJp-E' ' TT J„ LUp-E' 



TT ./o 



Qp Y^cou (-^ £;')di;' 



7^7 



Gb 



(30) 



where Cp is the integral from the high energy cutoff, f7p, to infinity and over the negative energies which give only 
weak energy dependence. 

In order to better understand the from of Sf°"* (p, E) near E ~ 3A, 5]^°"* {p, E) and S^°" (p, E) can be approxi- 



7^7 



7^7 



mated by simple analytic forms. From Figure 3 these can be approximate by S'^'J"* (p^ E) — ap{E —?iA)Q{E —ZA) with 



7^7 



7'7 



ap ~ 0.055. Xl^°" (p, E) increases more rapidly than this linear function by _B ~ 5A but this approximation is suffi- 



7'7 



cient to give a semi-quantitative account. I take S™" (p, E) — l3pO{Ep+3A — E)&{E—Ep)f{E—Ep), where f{E—Ep) 



7T7 



is a function approximately the form of the energy dependence of the collective mode contribution. Looking at the 



form of the collective mode contribution in Figure 3 it can be approximated hy f{E — Ep) ~ { 



_ /(E-£p)\2 



h3A-B 



Taking this form and requiring that the maximum value of S"^" (p, E) be 0.09A as in Figure 3 requires (3p ~ 0.02A 
for Ep > 3A. (3p grows from zero at _Ep = A and is independent of Ep beyond Ep = 3A. 
Putting in the forms of the two contributions to the S"^ (p, E) 



^'Jp,Ep)=.'^ 



^p {E' - 3A)dE' (3p 



3A 



Ep-E' 



S„+3A 



f{E')dE' 
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^^^^Q SA) + ^{Ep 
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(31) 



and 



E'^,^{p,up)^j:'^,^{p,Ep) + ^ 



TT 



2(3) 



3/2 fp_ 



Ll!„ — Ep 



A 



{ujp — 3A)ln 



V3 h^p^ 



3A| 
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3A 
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(32) 



From Figure 11 it can be seen that (5E' ^ (p, ujp) ~ 0.01 A so that tjp — £'p ~ S' ^ (p, E) giving 



7^7^ 



^\(p,ujp)o,j:'(p,Ep)il + 2{3) 
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(33) 



14 



Curve A:z: -H (p,E )/A vs- 



•2 




Curve B:(co — E )/z^ vs. cjl> /A 



^nergy(in units of ^} 



FIG. 11. The on-shell self-energy, , as a function of Ep and the shift in the quasiparticle spectrum, uip 
from Ep as a function of ujp for e^pp = 25.13A~^, uju = O.leF and A = 0.0273e_F. The position of the 
features in the quasiparticle spectrum are seen to coincide with the frequencies at which the strong-coupling 
features are seen in the density of states. 

Using (3p discussed above one finds from Equation (33) that the magnitude of E' ^ {p,ujp) is enhanced over that 
of E' I {p, Ep) by about 5%. Considering the simple form of the functions used this is consistent with the results 

shown in Figure 11. Taking Op and ap to be constants and allowing /3p to increase by about | between Ep — 2 A and 
for Ep — 3A, consistent the calculation of E"^ {p, E), the form of E' ^ (p, Ep) for Ep > 2.5A is seen to come from 

the combination of the ^{Ep — 3A)/n ' ''^ — - and the ^ ^"^^^ terms. The weak features in curve B, {ujp — Ep)/ IS. 
vs.wp/A, at Ldp = 3A ± jE' ^ (p, Ep) come from the term proportional to ap in Equation (34) where the term in 
the brackets goes through maxima in magnitude. However a better approximation for the momentum and energy 
dependence of E ^ (p, E) is required beyond the form used here in order to recover the Ep dependence of E' ^ (p, Ep) 
between Ep = 2 A and Ep — 2. 5 A without a detailed fitting of the Cp and Op functions. The (3p term coming from 
the collective contribution to E ^ (p, E) is responsible for the bump feature in E ^ (p, Ep) between 2.6A < E < 3A. 
It is clear that without the collective mode contribution the bump in the density of states below to = 3A would be 
absent and the step at uj — 3A would be more pronounced. 

VI. DISCUSSION 

As well as tunneling conductance and the density of states discussed here ARPES experiments can provide infor- 
mation on the single-quasiparticle spectral density once background has been taken into account. In the calculation 
discussed above the spectral density is a Lorentzian with no experimentally detectable non monotonic energy de- 
pendence away from the peak. Many-body effects give a non-trivial quasiparticle spectrum, ujp, leading to the 
strong-coupling effects discussed above but do not lead to any signatures in ARPES measurement such as the dip 
feature at w = 2A in the cuprate superconductors. A necessary condition that features away from the quasiparticle 
peak in the spectral density lead to observable features in ARPES is that the imaginary part of the self-energy should 
have both a strong energy dependence and be sufficiently large that there is weight away the peak. In the present case 
the linear dependence of E".f {p, Ep) on Ep and the small magnitude combine to ensure that the energy dependence 
of T,'\ (p, E) leads to no features in ARPES or tunnehng conductance. 



15 



This is in contrast to models used to describe ARPES and tunneling data in the cuprates |g|,^,|9J where the 
strong-coupling feature is at 2A in the density of states and the spectral density, as measured in ARPES, [0 and 
at 3 A in the SIS tunnehng conductance. In these models a phenomenological form for T,^f^{p,E) is used in which 
there is a rapid increase in the magnitude of E"^ {p, E) to a value '^ A in some momentum directions at E ~ 2A 
after which any dependence of S''^ {p, E) on E is weak. In these models E' ^ {p, E) is assumed to lead only to a 
renormalization of the gap. This model for E''^ (p, E) in the case of a d-wave superconductor is suggested by Golden- 
rule calculation of the decay of mean field quasiparticles in a tight-binding bandstructure P,0|. The results of the 
tight-binding calculations depend on both the chemical potential and the presence of next-nearest neighbor terms 
in the bandstructure. The form of E''^ (p, E) is consistent with that arising from low energy antiferromagnctic spin 
fluctuations l45|. The bandstructure and short-range correlations leading to low energy magnetic excitations which 
are responsible for the strong-coupling effects in the models of the cuprates are absent in the calculation reported 
here. However the present calculation points to the necessity of including long range Coulomb interaction in order to 
distinguish low dimensional from conventional superconductors and suggests that any model which is to quantitatively 
account for these features in the cuprates should incorporate the effects of long-range correlations as well as the short- 
range correlation associated with the magnetic properties of the cuprates. The extent to which a superconductor 
is low dimensional, where the effects discussed would be important,! is determined in the model described here by 
the weight in the effective quasiparticle interaction due to the Coulomb interaction at low frequencies. In layered 
superconductors Das Sarma et al. i29| , |30|l have shown that depending on interplanar coupling and on the direction 
of the wavevector the collective mode spectrum, cj^-, is at energies ^ A as |(7| ^ 0. There should still be substantial 
weight in the collective mode at frequencies < lOA in order for quasiparticle interactions to be enhanced above the 
3D case. This issue is presently being investigated p^ . 

Since the original identification of the dip feature in ARPES there have been suggestions that it is strongly associated 
with the mechanism for superconductivity in the cuprates. Arnold et al. [}46| identified the strong-coupling features 
seen in the data of Dessau et al. [16|. with a "bosonic peak" which they proceeded to analyze in terms of the a^F{ijj) 
function of Eliashberg theory. Analysis of more recent data have lead Shen and Schrieffer [|8[ to a model where the 
quasiparticles are strongly coupled to a collective mode. This collective mode contributes to the line shape in ARPES 
in the normal state and is strongly influenced by the superconducting transition. In their picture the dip feature 
in the line shape is the consequence of the strong quasiparticle coupling to the collective mode which provides the 
pairing mechanism for hiTc superconductivity. In the present case also the collective mode is strongly effected by the 
transition to superconducting ground state although it is not associated with the mechanism for superconductivity. 
It decays into pairs of superconducting quasiparticles once its energy reaches 2A, reemerges at higher energies and 
essentially follows the dispersion of the normal state plasmon. The contribution to off-shell self-energy, E (p, E), from 
scattering from the plasmon is very weak at the low energies, ~ A, for the 2D electron gas. It leads to a weak feature 
in the spectral density even at p = l.lp_F which corresponds to Ep ~ 8 A for A = 0.0273eF and would not be detected 
easily in experiment [E^. 

VII. CONCLUSIONS 

In this paper I have shown how quantum fluctuations in 2D charged Fermi systems are enhanced over those in 3D. 
This result holds for both the normal and superconducting states. However in the superconducting state they clearly 
lead to strong-coupling features in experimental quantities which are a probe of the superconducting state. 

In the model discussed in this paper the quasiparticle interaction responsible for the feature at w = 3A in SIN for 
s-wave superconductor and the corresponding feature at w = 4A in SIS is the Coulomb interaction and not the pairing 
interaction giving the s-wave superconducting ground state. In the absence of the Coulomb interaction the pairing 
interaction would produce the same feature but would not account for the difference in magnitude of these effects in 
conventional 3D superconductors and the quasi-2D cuprate superconductors. The anomalously large magnitude of 
the strong coupling features seen in the cuprates suggest that the electronic properties of these materials are quasi-2D. 
This is a basic assumption which most workers in the area have employed for a long time. However the long range 
part of the Coulomb interaction has been ignored. The present calculation suggests that, although the interactions 
which provide the mechanism for superconductivity in the cuprates may be short correlations, long range Coulomb 
correlations should also be included. In particular the pairing mechanism undoubtedly contributes to the magnitude of 
the feature at w = 2A in ARPES data and to the strong-coupling features seen in tunneling experiments at eV — 3A 
but they are unlikely to solely account for its magnitude. 
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